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Abstract 

In a recent experiment, Lai et al. [Phys. Rev. B 75 (2007) 033314] studied the apparent metal-insulator transition (MIT) of a Si 
quantum well structure. Tuning the charge carrier concentration n, they measured the conductivity cr(T, ri) for a very dense set 
of n values. They observed linear T dependences of cr around the Fermi temperature and found that the corresponding T — > 
extrapolation cr () (n) exhibits a sharp bend just at the MIT. Reconsidering the data by Lai et at, it is shown here that this sharp 
bend is related to a peculiarity of <t(T = const., ri), which is clearly detectable in the whole T range up to 4 K, the highest 
measuring temperature in that work. It may indicate a sharp continuous phase transition between the regions of apparent metallic 
and activated conduction to be present at finite temperature. This interpretation is confirmed by a scaling analysis without fit, which 
illuminates similarities to previous experiments and provides understanding of the shape of the peculiarity and of sharp peaks found 
in dlog 10 cr/dn (ri). Simultaneously, the scaling analysis uncovers a strange feature of the apparent metallic state. 

Key words: Metal-insulator transition, localisation, scaling, apparent metallic phase 
PACS: 71.30.+h, 73.20. Fz, 73.40. Qv, 73.63. Hs 
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1. Introduction 

Is electronic transport in 2-dimensional (2d) disordered sys- 
tems exclusively nonmetallic, so that the resistivity p(T) always 
diverges as the temperature T tends to 0, or can it also have 
metallic character? This fundamental problem has been con- 
troversially debated for three decades: The existence of a cor- 
responding metal-insulator transition (MIT) was denied by the 
localisation theory of Abrahams et al. [1], which, however, ne- 
glects electron-electron interaction. Thus, it came as a surprise 
when Kravchenko et al., who had varied the charge carrier con- 
centration n in high mobility MOSFET samples, first reported 
a strong decrease of p down to 20 mK [2]. They considered the 
conduction in the respective (T, ri) region as metallic, an inter- 
pretation which was questioned in particular by Altshuler and 
Maslov [J]. The nature of the apparent metallic state has re- 
mained puzzling up to now, for recent reviews see Refs. IHIH]. 

Recently, Lai and coworkers studied the transport properties 
close to the apparent MIT of an n-type Si quantum well con- 
fined in a Sio.75Geo.25/Si/Sio.75Geo.25 heterostructure [6]. The 
authors observed the conductivity cr to have nearly linear T de- 
pendences around the Fermi temperature Tp varying between 
2 and 2.5 K. They showed that <x (n), the T — conductiv- 
ity obtained by linear extrapolation from this T region, exhibits 
two regimes of different slope. Particularly important, there is 
a sharp bend at the transition between the two regimes. It co- 
incides with the n value, n , where dcr/dT changes its sign as 
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T — > 0. The authors interpret this finding as an indication of 
the existence of two different phases. Thus, the apparent MIT 
should be related to a real phase transition at T — 0. This tran- 
sition is commonly expected to be smoothed out at finite T. 

However, such a smoothing is not a must for a quantum phase 
transition as is obvious from two examples: The Ising ferro- 
magnet LiHoF4 in transverse magnetic field undergoes a quan- 
tum phase transition at T = 0, where the quantum critical point 
is the beginning of a line of continuous phase transitions at finite 
T Q711 . The phase diagram of the quasi-two-dimensional organic 
conductor k-(BEDT-TTF)2X exhibits a line of first-order transi- 
tions at finite T separating an antiferromagnetic-insulator phase 
from an unconventional-superconductor phase at low T [81]. 

From this perspective, the additional observation by Lai et 
al., that the slope of cr(T) for T ~ Tp is almost constant 
close to the MIT, n = n c , provokes two interesting questions: 
(i) Provided the extrapolated cro(ri) has a knee indeed, and 
the slope used in extrapolation is roughly constant, should not 
cr(T = const., ri) exhibit a knee at « c also for finite Tl (ii) If 
yes, is the existence of this peculiarity restricted to the region 
of linear cr(T) around Tp, or is it a more general phenomenon? 

Starting from these questions, the present work points to sev- 
eral indications for the MIT in 2d systems being connected with 
a line of continuous phase transitions at finite T. It substan- 
tially extends the detailed arguments given in Ref. [9]: More 
data are taken into consideration in the scaling region, where 
cr(T,ri) = cr((« - n c )/T l IP), consequences of this scaling for 
the physical nature of the apparent metallic state are discussed, 
and possible corresponding similarities with the transport in an 
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Figure 1 : Charge carrier concentration dependence of the conductivity for the 
Si quantum well structure studied in Ref . [ 6] . Data were obtained from Fig. 1(b) 
of that work. The arrow marks the critical concentration n c of the apparent MIT, 
see text. The dashed line is a fit by the piecewise linear function introduced in 
the text. The full lines serve as guide to the eyes. 



AlAs quantum well [ 10] are illuminated. 



2. Phenomenological reconsideration of the experiment by 
Lai and coworkers 

To find out to which extent and in which form the sharp bend 
of <r (n) persists when T is finite, I digitised the high-resolution 
preprint version of Fig. 1(b) of Ref. [|6[]. Fig. 1 depicts the ob- 
tained cr(T = const., n) data for five T values. First of all, this 
graph shows that all cr(T = const., n) curves share a peculiar 
feature, namely an "indentation" at n « 0.32 • 10 11 cm 2 . It oc- 
curs close to the concentration value where dcr/dT changes 
its sign at the lowest T considered in Ref. [6]. This value, 
n c - 0.322 ■ 10 11 cirT 2 , is the critical concentration of the ap- 
parent MIT. It is marked by arrows in Figs. 1, 2, and 4. 

To answer the question whether or not the sharp bend ob- 
served by Lai et al. at <r («) persists at finite T, fits by the piece- 
wise linear function f p \f(n) = a + bn + c(n - n^)6(n - «k) with 
the adjustable parameters a, b, c, and n\ were tried. Here, «k 
denotes the knee position, and 9 stands for the Heaviside step 
function. For comparison, the data sets were approximated by 
polynomials of third order, which have the same number of ad- 
justable parameters. Additionally to the data shown in Fig. 1, 
<t(0.47K, n) and cr(0.87K, n) were included in this analysis. 

Table 1, left part, presents the optimum knee positions 
and the related values of the least square sum^- 2 , where weight 
1 was ascribed to all points. It is remarkable that the same ny 
value, 0.318 ■ 10 11 cm" 2 , results for all temperatures but 0.47 
and 0.35 K. This «k value is only very slightly lower than the 
abscissa of the knee for the T — * extrapolation cr () (n) in Fig. 2 
of Ref. |0]. Moreover, note that, for the four highest T values, 
piecewise linear functions approximate the experimental data 
clearly better than the polynomials. 

The advantage of piecewise linear functions over polynomi- 
als of third order diminishes with decreasing T. This trend 



Table 1: Results of phenomenological fits of cr(T = const., n) data from the Si 
quantum well study Ref. 1 6] . Piecewise linear functions (plf) and polynomials 
of third order (pto) are compared. Values of knee positions % and least square 
siims^ 3 are given in units of 10" cm~ 2 and e 4 /h 2 , respectively. The left group 
of columns results from consideration of the complete data sets, the right group 
is obtained from the data points fulfilling 0.5 <t(T, n c ) < cr(T, n) < 2.0 <r(T, n c ). 
N denotes the number of data points taken into account. 
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is not surprising because, in the region of activated conduc- 
tion, <t(T = const., «) vanishes in some exponential way with 
decreasing n. Hence, when considering a fixed n range and 
lowering T, the cr region broadens, and the basic nonlinear- 
ity gets increasing weight compared to the "indentation" at 
n w 0.32 • 10 11 cnT 2 . Thus, a more meaningful comparison 
of both approximations is obtained restricting all fits to compa- 
rable cr intervals, here to {0.5 <x(T, n c ), 2.0 cr(T, n c )}. Under this 
condition, for each of the T values considered, piecewise linear 
functions approximate the data clearly better than polynomials 
of third order, see Tab. 1, right part. Moreover, now all fits yield 
the same value, 0.318 • 10 11 cirT 2 . 

Thus, for all T values, some characteristic change occurs 
within a very small n region close to the knee of the piecewise 
linear fit. In this sense, the sharp bend observed by Lai et al. 
at the T — » extrapolation <T (n) turns out to be not smoothed 
for finite T. In particular, it seems to exist even for T <k Tp, 
where <x(T, n = const.) is nonlinear and considerably deviates 
from the linear extrapolation estimate from the vicinity of Tp. 
Because « « c for all T, it is likely that the bends originate 
from sharp continuous phase transitions separating the regions 
of apparent metallic and activated conduction at finite T. 

However, there is a tiny systematic difference between the 
values of ny and n c . It presumably arises from the approxima- 
tion by piecewise linear functions being an oversimplification. 
This aspect is considered in the next section. 



3. Scaling analysis 

Provided the sharp bend is indeed related to a phase transi- 
tion, how may it arise? What fine structure could it have? On 
a macroscopic level, o~(T, n) can be described in the following 
way. Suppose, on the insulating side of the apparent MIT, the T 
dependence of cr scales as observed by Rravchenko et al. Ill 1 1 
at the resistivity of MOSFETs, 



cr(T, n) = cr(f) with t = T/T (n) 
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Figure 2: Scaling check without adjustment of parameters considering the four 
<t(T = const, n) data sets with T < IK which are analysed in Tab. 1. 



For further references on scaling of T dependences close to the 
MIT in various 2d and 3d systems, see |9{]- Assume that this 
scaling holds up to n c where 7q(h) — > 0. Suppose furthermore 



7o(n) = A\6nf with 5n = n 



(2) 



where A is a constant, which might, however, be sample depen- 
dent. Thus, cr should depend only on the quotient T/\6nf or 
rather on Sn/T 1 ^, compare Fig. 10 of Ref. 11 ill . 

The value of the critical exponent /3 should be universal. It 
was determined by Kravchenko et al. ifTHl from their MOSFET 
experiment, /? = 1.6 ± 0.1. The value n c = 0.322 ■ 10 11 cm -2 
is known from the sign change of d cr/d T as T — > 0. Thus, 
without adjusting any parameter, a scaling check of cr(T, n) can 
be performed plotting <x as function of Sn/T 1 ^ for different T, 
This is done in Fig. 2 considering the <x(T = const, n) data sets 
for T = 0.35, 0.47, 0.64, and 0.87 K. Data sets for higher T 
are not taken into account since <t(T, n c ) becomes T dependent 
above roughly 1 K, presumably because some additional mech- 
anism becomes relevant there, see Fig. 1 of Ref. Jg]. 

Note that such a scaling check ascribes higher weight to the 
points more distant from n c , while the analysis in the previous 
section is particular sensitive with respect to the data very close 
to n . Therefore, here, the influence of sample inhomogeneities 
is reduced in comparison to the fits in Tab. 1 . 

Figure 2 shows that, for n < n c , all four o-(6n/T 625 ) curves 
nicely fall together although no parameter has been adjusted. 
It is unlikely that the agreement of the curves arises only by 
chance, so that Fig. 2 can be considered as support of the scaling 
hypothesis. Hence, since To(n) — * as n — * n c , this graph 
indirectly indicates phase transitions at finite T. 

Also for n > n c , that means within the apparent metallic 
phase, the four curves in Fig. 2 collapse. This is surprising: In 
case of conventional metallic conduction, as T vanishes, cr(T, n) 
would tend to a finite and n dependent value, <x(0, n), which in- 
creases monotonously with n. Vanishing T corresponds to di- 
verging 6n/T - 625 , so that curves, which were drawn for varying 
T and fixed 6n in Fig. 2, would split, in contradiction to the ob- 
served scaling. Thus, scaling for n > n c cannot be understood 
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Figure 3: Scaling check for AlAs quantum well data from Fig. 2 of Ref. 1 10]: 
(a) high-mobility and (b) low-mobility directions. Here, n c = 0.70 ■ 10 11 cm" 2 
is presumed, see text. 



in terms of conventional metallic behaviour. Superconductivity 
at T — might be an alternative, compare the discussion in Ref. 
1 1 ill and the study of ultrathin Bi films by Liu et al. 11211 . 

However, for n > n c , a lot of publications have reported sat- 
uration of cr as T — > 0. Thus several questions arise: Might 
the scaling observed here be restricted to the considered ma- 
terial or to a tight vicinity of n c ? What about the influence of 
inhomogeneities, sample size or thermal decoupling? 

To get a preliminary answer to the first question, reconsider 
now an experimental study of an AlAs quantum well by Pa- 
padakis and Shayegan [10]. Figure 3 shows scaling checks 
for data obtained by digitising Fig. 2 of Ref. [10]. To cor- 
rect for the T dependence of cr at n c , the quotient q(T, n) = 
<t(T ', n) I <t(T ', n c ) is considered here instead of <r(T,ri) itself. 
This approach turned out to be effective for the Si quantum well 
data by Lai et al. in Ref. |@]. 

The scaling check given in Fig. 3 has a hidden parameter: 
Contrary to the situation in Fig. 2, the n c value has some uncer- 
tainty. According to Fig. 2 of Ref. [10], the « c values defined 
by dcr/dr = as T -» are roughly 0.65 and 0.70 ■ 10 11 cm" 2 
for the high- and low-mobility directions, respectively, but n c 
should be uniquely defined. (A possible origin of the direction 
dependence of the mobility is discussed in detail in Ref. llOfl.) 
In comparing plots for several values of n c , the best data col- 
lapse was observed for n c = 0.70 • 10 11 cm" 2 . 

Figure 3 (broader n range) resembles Fig. 2 to a large ex- 
tent. Thus, although the random dispersion of the data points is 
here clearly higher than in Fig. 2, the proposed scaling approach 
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Figure 4: dlog 10 cr/d« versus n for the Si quantum well curves displayed in Fig. 
1. Concerning details of the calculation of the derivative values see Ref. (^|. 

might be applicable to the AlAs data from Ref. fioll . Moreover, 
also in Fig. 3, the slope changes significantly at n — n c . 

Turn now again to the analysis of the Si quantum well data 
by Lai et al. The "indentation" around 5n — in Fig. 2 can be 
understood in the following way. Since 7o(n) — > and t —* oo 
as n — > n c for arbitrary constant T, the ansatz 

cr(*) = ov(l-Bf v ) (3) 

is suggested. Here B is a dimensionless constant and v a posi- 
tive exponent. Hence, 

o-(T,n) = cr c -(l-C\6nf v ) with C = B(A/T) V . (4) 

Only a rough guess of the value of v can be given at present. 
In case of hopping in the Coulomb gap for small t, v would be 
1/2. However, according to Ref. [11], v may be expected to 
be smaller for large t. In both cases, the product Bv would be 
clearly smaller than 1 so that cr(T = const., n) should have a 
root-like peculiarity at n c . 

This hypothetical peculiarity implies a divergence of dcr/dn, 
and thus also of dlog 10 cr/dn. Already at finite temperature, it 
indicates the transition between activated and apparent metallic 
conduction. As a consequence, it permits the precise identifica- 
tion of the MIT. This approach resembles the strategy used in 
studying ultrathin metal films in Ref. 11311 . where the crossover 
between logarithmic and exponential temperature dependence 
was related to a change in the thickness dependence of the 
thickness derivative of the conductance. 

In reality, the expected divergence would be smoothed to a 
sharp maximum. Fig. 3 shows that, indeed, d log 10 cr/d n ex- 
hibits such a sharp peak independent of T - note that a logarith- 
mic scale is used to display the derivative values. It is striking, 
although a natural consequence of Eq. (4), that the peak is al- 
ways located just at n c , another indication of a phase transition 
occurring at finite T. Of course, the question arises whether 
or not this sharp maximum is specific to the Si quantum well 
studied in Ref. [60. The comparison of Fig. 2 with Fig. 3 sug- 
gests that also the AlAs quantum well studied in Ref. 1 10] may 
exhibit such a feature. 



4. Challenges 

As a whole, the presented indications give strong support to 
the hypothesis that the MIT at T — is connected with a line of 
sharp phase transitions at finite T. However, this conclusion has 
to be confirmed by further experiments. Such studies should 
primarily focus on enhancing precision, in particular concern- 
ing homogeneity of the samples, rather than on reducing the 
lowest accessible T for the following reasons: In plots anal- 
ogous to Fig. 2, the effect of inhomogeneities increases with 
decreasing T. For arbitrary T, inhomogeneities smooth the in- 
dentations shown in Figs. 1 to 3 and the peaks in Fig. 4. 

In case, the hypothetical line of phase transitions at finite T 
exists indeed, three questions arise: Does it terminate at a cer- 
tain T value? Under which conditions does scaling hold in the 
region of apparent metallic conduction, is there a specific (T, n) 
region? What is the nature of the related unconventional phase? 

I am indebted to T. Vojta for an important literature hint. 
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